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O 

. We present a decomposition of the hadronic tensor for a general polarized Drell-Yan process 

AB — ¥ £ + £~ X in terms of helicity structure functions using the helicity axes of the dilepton rest 
frame as a basis. Next, we consider a QCD parton model and in the framework of a generalized QCD 
factorization scheme, which applies when the dilepton invariant mass M is much larger than the 
transverse component of the photon momentum, qr, in the hadronic center of mass frame. In this 
approximation we compute the angular distribution of the unpolarized Drell Yan cross section and 
of the Sivers effect, in single polarized Drell-Yan scattering, by taking into account the transverse 
motion of partons inside the initial hadrons, kj_. Interesting and simple results are found in the 
(-H ■ kinematical region k± ~ qr -C M, to first order in a qr/M expansion. Finally, explicit analytical 

expressions, convenient for phenomenological studies, are obtained assuming a factorized Gaussian 
dependence on intrinsic momenta for the unpolarized and Sivers distribution functions, in analogy 
to those obtained for semi-inclusive deep inelastic scattering. 
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Lepton pair production in hadronic collisions, AB — > £ + £ X, is one of the most useful tools to study 
hadronic interactions and is presently the focus of considerable interest from both experimental and theoret- 
ical points of view. Drell-Yan (DY) scattering processes in facts, due to their purely leptonic final state, can 
be considered among the cleanest processes in order to study the hadron structure, including the puzzling 
question of the nucleon spin structure. Moreover, at Born level, they are purely electromagnetic processes, 
involving one elementary scattering channel {qq annihilation): therefore they play, in hadronic interactions, 
the same, crucial role of DIS processes in electromagnetic interactions. 

Here we present a decomposition of the hadronic tensor for a general polarized DY process in terms of 
"helicity" structure functions. Similar decompositions were presented in Ref. [l| for the unpolarized case 
and in Ref. Q for the single polarized case. In Ref. Q the full expression for the general polarized DY 
cross section was obtained by using "invariant" structure functions, i.e. using the spin vectors, the hadronic 
momenta and the dilepton momentum as a basis. This method is rather technical and lengthy. Instead, we 
will show that the use of "helicity" structure functions, inspired by the approach presented in Ref. |l|, is 
considerably simpler, elegant and easy to check. Moreover the helicity structure functions can be directly 
related to hadronic degrees of freedom, i.e. they can be given a partonic interpretation. Invariant structure 
functions and helicity structure functions are obviously related to each other, as explained in details in 
Ref. |l|. We present some relevant results obtained in a non-collinear extension of the QCD parton model 
framework, where we take into account the transverse motion of the constituents: the kinematical configura- 
tion of the process becomes much richer, allowing for more degrees of freedom and new possible correlations 
between hadron and parton spins and momenta. Thus, beside the three usual collinear distribution functions 
(unpolarized, helicity and transversity distribution functions), five more Transverse Momentum Dependent 
distribution functions (TMDs) are necessary to describe the nucleon structure. Among them we find the 
Sivers function 0, [f| , which gives the probability to find an unpolarized parton inside a transversally polar- 
ized nucleon, and the so called Boer-Mulders function describing a transversally polarized parton inside an 
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unpolarized nucleon. The TMDs usually exhibit highly non-trivial QCD properties (for instance, some of 
them are naively T-odd). Their theoretical definition and their QCD evolution are still largely unknown. In 
spite of that, there are many experimental evidences on the existence of these phenomena, provided by sev- 
eral experimental collaborations in many facilities in the world: HERMES at DESY, COMPASS at CERN, 
RHIC at BNL, BELLE at KEK, and JLAB (for a theoretical and experimental comprehensive review, see 
Ref. [6| and references therein). 

The inadequacy of the simple, collinear QCD parton model emerging from the latest experimental data 
involving polarized hadronic reactions motivates our efforts towards the study of a more refined picture, 
in which the "transverse" , three-dimensional structure of hadrons relevantly comes into play, offering a 
fundamental ground for the investigation of the detailed form of intrinsic transverse momentum dependent 
(TMD) parton distribution functions (PDFs), which embed all the complexity of the above-mentioned 3-D 
structure of hadrons. 

While in inclusive DIS the parton intrinsic transverse momentum of partons is invisible in any observable, 
as it is totally reabsorbed in the integration, this does not happen in processes like Drell-Yan or Semi- 
Inclusive DIS (SIDIS), where it becomes manifest in certain observables, like the angular distributions of 
the lepton pairs (i.e. the cross sections differential with respect to qx in DY) or the angular distributions 
of the final detected hadron (i.e. the cross sections differential with respect to Pt, the hadronic transverse 
momentum, in SIDIS). 

In Sect. U we start by giving a brief review of the general decomposition of the unpolarized DY cross 
section in terms of helicity structure functions [l|, which we prefer to the invariant structure functions as 
they naturally show the hadronic structure properties in terms of their elementary constituents. While 
performing the few necessary steps to obtain the hadronic tensor decomposition, we set up a convenient 
methodology which will naturally lead our way through the much more complex polarized case. 

In Sect. |TT] we define a covariant hadronic spin vector, and decompose the single polarized term of the 
hadronic tensor. We conclude with some interesting comments about the angular dependence of the single 
polarized DY cross section in two of the most used reference frames, the Collins-Soper (CS) dilepton rest 
frame 0] and the hadronic center of mass (cm.) frame. Similarly, in Sect. lIIll we work out the decomposition 
of the doubly polarized term of the hadronic tensor and obtain the corresponding angular distributions. 

In Sect. [TV] we work in the QCD parton model approximation, where factorization allows us to write the 
DY hadronic tensor as the incoherent sum of the partonic tensors of the corresponding parton annihilation 
processes. For the time being, we consider only the contributions to the DY cross section generated by 
unpolarized qq — > £ + £~ elementary cross section, i.e. the terms proportional to the convolution of two 
unpolarized TMD PDFs and to the Sivers effect. The remaining polarized terms will be considered elsewhere. 
The full angular distribution of such terms is given in the CS frame, to first order in a qx/M expansion 
(where qx is the photon transverse momentum as measured in the hadronic cm. frame). Finally, the 
phenomcnological implications of our results are discussed by considering a simple Gaussian model for the 
k± distribution of the TMDs. 

Our conclusions are summarized in Sect. |Vj 



I. BRIEF REVIEW OF THE UNPOLARIZED CASE 



The general expression for the angular distribution of the polarized and unpolarized Drell-Yan process 
can be predicted by means of a kinematical analysis. The cross section for the process AB — > £ + £~X can 
be expressed in terms of the hadronic and leptonic tensors, W^ u and L^ u respectively, [l|, @] as: 

da a2 ^W" (1) 



d 4 qdn 2(2tt 4 )sM 4 ' 



where q = (<7o> <Zt> <?£) is the dilepton momentum in the hadronic cm. frame, while q 2 — M 2 is the dilepton 
mass squared. The unpolarized, symmetric leptonic tensor is well known and can easily be written in term 
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of the lepton momenta, I3 and I4: 

L" v = 2%ll + 2l%l% - M 2 g^ (2) 

The hadronic tensor is a non-perturbative object describing the dynamic of the interaction between the two 
colliding hadrons. However, it is not totally unknown as it must fulfill fundamental constraints like: 

q l iW lv = q v W^ u = current conservation; (3) 

w ^u = hermiticity; (4) 

-pW^y- 1 = W^ v parity invariance. (5) 

where V is the parity operator. Moreover, the hadronic tensor must be a function of the physical invariants of 
the process. Consequently, it is possible to decompose the hadronic tensor by means of structure functions, in 
complete analogy to what is usually done in deep inelastic processes. In particular here we want decompose 
the hadronic tensor by means of the so-called helicity structure functions, i.e. projecting the hadronic tensor 
on the components of the virtual photon polarization vector [JJ . 

For the unpolarized Drell-Yan process this leads to the well known decomposition in terms of the unpo- 
larized helicity structure functions Wt, Wl, Wa, Waa 

W^ = (-g^ + (W T + Waa) ~ 2x»x u Waa 

+ z^z u {W L -W T - Waa) ~ [x^z v + z^x")W A (6) 

where Wt and Wl are the structure functions for transversely and longitudinally polarized virtual photon, 
while Wa and Waa correspond to single and double spin flip respectively; x^, z^ (and y M ) are the dilepton 
rest frame axes. They are four- vectors defined in such a way that they are orthogonal to each other and to 
the four- vector q^, and transform in unitary three- vectors in the dilepton rest frame. For convenience the 
axes are normalized to —1. An invariant expression of x^ , y^, z^ can be found in Ref. Q . Notice that 
W£ = ~(2W T + W L ), having followed the convention g 00 = I, g u = -1 (i = 1,2,3) while Ref. Q follows 
the opposite convention. Contracting the hadronic tensor with the leptonic tensor, Eqs. (f2l IA4|) . one obtains 
the most general structure for the unpolarized Drell-Yan angular distribution: 



da unp a 2 



<m>d 4 q 2(2ir) 4 M 2 s 

+ W A sin 29' cos 4>' + Waa sin 2 9' cos 20' 



W T (1 + cos 2 9') + W L (l - cos 2 9') 

(7) 



where 9' and are respectively the lepton polar and azimuthal angles in the dilepton rest frame defined 
by the chosen axes x M , y^, z 11 . A similar expression for the single polarized case can be found in Ref. Q 
while Ref. Q shows the most general case where both hadrons can be polarized. Here we present a different, 
elegant derivation using the properties of the hadronic tensor, Eqs. (|3|4l5j) . and a decomposition of this 
tensor in terms of helicity structure functions projected onto the helicity axes of the dilepton rest frame, in 
complete analogy to the unpolarized case outlined above. In order to explore the polarized case it is better 
to use the more flexible notation of Ref. For example, we rewrite Eq. ([6]) in this way: 



W^ = (-<T + iJLJ (Fhu + F™ 2 *') - 2x^F^' 

+z»z» (F 2 V - Fbu - F^') - (x»z» + Z W) F™*' (8) 



where UU denotes that both hadrons are unpolarized and F^y = Wt, F^ v = Wl, = Wa and 

^,ras20 _ Waa- This notation has the advantage of reminding the azimuthal dependence of each helicity 
structure function. 
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II. SINGLE POLARIZED DRELL-YAN 



To deal with the polarized Drell Yan processes, we need to properly take into account the polarization 
degrees of freedom. For this reason, it is convenient to define the following four-vectors: 

&h = (^o/n^r/i)Sifi) (9) 

S£ h = (0,S Th ,0) (10) 

S£ h = (S oh ,0,S Lh ) (11) 

= (z-S h )z» (12) 

S lh = Sjl-(^-+z^)s hl/ , (13) 

where S% = {Soh, Sxh, Slh) is the covariant spin vector of hadron h — A, B in the hadronic cm. frame. It 
can be easily seen that S^ h q^ = 0, S^q^ — 0. This property is very important because any decomposition of 
the hadronic tensor by means of the helicity axes x^, y^, z^ and the S^ h and S± h vectors automatically fulfills 
the requirement of Eq. (|3]). Notice also that Sj_ h and Sjfu are purely transverse (0, S' ±h ,0) or longitudinal 
(0,0,0,5^) (three- vectors) in the dileptonic rest frame. 
In general, the hadronic polarized tensor can be written as a sum of four terms: 

W^(S A , S B ) = WZ P + + Wg + W^ aSb , (14) 

where W£% p is the unpolarized hadronic tensor of Eqs. ((U or (j8]); Wg V A and Wg V B are the single polarized 
terms of the hadronic tensor while Wg" s is the double polarized term. 

Applying current conservation, hermiticity and parity conservation, Eqs. ([3][5]), the single polarized 
hadronic tensor can be decomposed in this way: 

WZ = - {S±A ■ V) (Flu + F?* 2 *') {-gT + 

+ 2(S ±A -V)(l 
+ (S± A -y)(^Fi l 
+ (S ±A -y)F c ^'(x^z v + z^x v ) 
- F^'iS^y" + y»Sl A ) - F^'(S^ y " + y^SjU) 
+ {S iA -z)F^'(x»y» +y»x») 

+ (SxA-x)FZ*'(z» y » + y»z»)- (15) 

In the first four lines we recognize the the same patterns found for the unpolarized cross section of Eq. ([5]). 
while in the last three lines we find genuinely new structures, proportional to the spin vector components. 
For the F structure functions we have adopted a notation similar to that used in Eq. ((8]), although their 
explicit connection to the stated angular dependence will only become known once the contraction with the 
leptonic tensor has been performed. 

The general, single polarized DY cross section can conveniently be separated in two parts 

da{S A ) _ da un P da s * 

dn'd 4 q dn'd 4 q d£l'd 4 q ' [ ' 

where ^, d I q is the unpolarized cross section of Eq. ©. Contracting the single polarized hadronic tensor 
of Eq. (fl"5l) with the leptonic tensor of Eq. (|A4|) . according to Eq. (QJ, we obtain the polarized term of the 
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F 2 , j-iCos2<^' 



5 



above cross section 



da SA 



dn'd 4 q 2(2tt) 4 M 2 s 



xl S' ±A sin0' SA 



+SW cos <t>' SA 



Fluil + cos 2 0') + Flu(l ~ cos 2 9') 

sin 29' cos 0' + F™* 20 ' sin 2 9' cos 20' 
Fj 1 " *' sin 20' sin 0' + F% 2 ^' sin 2 0' sin 2<ff 



+S' 



(A 



F™ sin 26' sin 0' + F™ 20 ' sin 2 0' sin 20 



(17) 



Expressions similar to Eqs. f| 15[) and (|T7|) can easily be written for polarized hadron B. Here the primes 
denote that the angles are in the leptonic cm. frame, analogously the polarization directions £ and _L in the 
F structure functions are referred to the dilepton rest frame, and 



(s: 



±A 
I 



cos< 



>J S A ) 

>'sa) 



- (Sj_ A ■ x) 

- {S±a ■ y) 



(18) 
(19) 



Notice that there is a straightforward correspondence between the polarization directions £ and _L in the 
leptonic rest frame and the longitudinal and transverse polarization directions, L and T, in the hadronic 
cm. frame. In particular 



sin ( 



whereas 



'S A . 



= -SIaVh = - S taVv = S ta sin( 



(20) 



5 



u 4 c 



_ CM 

M — "A^C 

— _ CM _ 
"At 



— CM , 

— d TA „ 



CM 



CM 



LA M ' 



and consequently 



(S^COS^gJ ^ S T cos(<j)s A ~ 7 ) 



(21) 



(22) 
(23) 



where 7 and 



S A are the dilepton and the proton A polarization azimuthal angles in the hadronic cm. 
frame. Notice that Eq. (|20l) holds because by definition the y^ axis is purely transverse in the hadronic cm. 
frame, while and are not. Therefore a transverse spin vector in the dilepton rest frame can acquire 
a longitudinal component in the hadronic cm. frame (and, vice-versa, a longitudinal spin vector in the 
dilepton rest frame can acquire a transverse component in the hadronic cm. frame). This was also pointed 
out in Ref. Similar expressions can be found in a slightly different notation in Ref. 0. 

For practical purposes, it is more useful to rewrite the polarized cross section of Eq. (fl7|) in terms of 
polarizations which are purely transverse and longitudinal in the hadronic cm. frame. This can be done 
inserting Eqs. (|20M21j) into the general expression of the hadronic polarized tensor Wg" of Eq. (|15|) and 
rearranging the structure functions. Here, we show how this can be done adopting the Collins Soper (CS) 
frame [7|, in which the kinematics is simple and transparent allowing to perform some useful approximations. 
Using the definitions of the CS axes x^, y^, as given in Appendix IB1 we obtain 



U TA^ ~ M 

s taVv = -Sta sin( 



Sta cos(4>s a 
Sa ~ 4h) 
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Sta z p- 



^ 7 
LA t 1 







qr (qo - qt) 



M sjM 2 + 



Sla 



(go - Ql) 
VM 2 + Qt 



Sla ■ 



(24) 



where qx and ql are, respectively, the transverse and longitudinal components of the virtual photon mo- 
mentum in the hadronic cm. frame. With these in mind, we can write the single polarized Drell Yan cross 
section as 

da SA a 2 



dWd 4 q 



2(2ir) 4 M 2 s 

x J St a sin 4> Sa \f tu (1 + cos 2 $')+F* u (i- cos 2 9') 



-StaFI 



TV 



sin 29' sin(0s A — 4> 1 — </>') 



+Sta-t T u sin 20 sm( 



+S T aF t Z { " S 

„sin(cj> s 
+JTA^TU 

+Sla 



h + 20') 



^7 +20') . 2 a' ■ Ix 

sin sm(05 A 



Flu sin 29' sin < 



F^ 20 'sin 2 0' sin 20' 



(25) 



having rearranged the i* 1 structure functions in the following way 

i(<t>s A -<t>-,T<t>') 



TU 



F 



, VM 2 + q T Train4 
M ^ 



F 



TV 



b LV 



^sin 2<^ 
LV 



j-icos 20 _i_ 
(90 - gL) 

(go - qi) 

V M2 + Qt 



j~ gj- Tjisin 20' 

M 



±(7 



gT j-isin <l 
~M U 



qT l? sin20' 



cv 



—F 

M 



v 



-F 



sin 2(f>' 



(26) 



It is interesting to notice that, in the Collins-Soper frame, in the limit qT <C M, one has (S'j_ a cos4>' Sa ) ~ 
cos((/>s A — 7 ), i.e. ~ Sta, and S^a — Sla- Therefore, in this limit, Eq. ([17]) holds with _L= T and 
£ = L, recovering the results of Ref. [3j], with </> 7 = 0. 



III. DOUBLE POLARIZED DRELL- YAN 



For the doubly polarized Drell Yan process we proceed in a similar way. For convenience, we decompose 
the double polarized hadronic tensor Wg^ s in four parts, according to the longitudinal and/or transverse 
components of the spin vectors Sa and Sb which they involve: 



w&\ =w%r q , +w%r „ +w&"„, « 

OAOB °lA a tB ^IJ^IB ^U^IB a ijl a (£ 



(27) 



By requiring current conservation, hermiticity and parity conservation, Eqs. ([3][5]), we obtain the decompo- 
sition of each term as a function of F helicity structure functions: 
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W»l s , B = (5 M • *) {SiB ■ z) (F}t + F™**') {-sT + 

- 2 [Sia ■ z) (SiB ■ z) (F^' + Ff* 2 *') x»x v 
+ (5m • z) {S m ■ z) (f* ~ F\ t - F^') z»z» 

- [S lA -z){S iB ■z)F™ s *'{x»z 1 ' + z»x"), (28) 



W S' ±A S> ±B = K S ±A ■ x ) (S±b ■ x) (S ±A ■ y) (S XB ■ y)} + F£**') (-<T + ^f) 

- 2 [(S ±A ■ x) (S ±B ■ x) - (S ±A ■ y) (S ±B ■ y)] (f^ a *' + x»x v 

+ [(S ±A ■ x) (S±g ■ x) - (S ±A ■ y) (S ±B ■ y)} (FL. - Fl ± - F ™ 2 *') z »z v 

- {(S ±A ■ x) (S±g ■ x) - (S ±A ■ y) (S ±B ■ y)} F c °f{x»z» + z»x v ) 

+ [(S ±A ■ x) (S ±B ■ x) + (S ±A ■ y) (S ±B ■ y)] + F cos2 ^') (-g^ + 

- 2 [(S ±A ■ x) (S ±B ■ x) + (S ±A ■ y) (S ±B ■ y)} (f™ 2 *' + F^') x»x v 



+ 


_A 


l) (5lB 


a;) 


+ {Sa_ a 


y) (S±b 


y)} (f! ± - f ±± 




- [(Sj 


_A 


l) (5_|_B 


x) 


+ (S ±A 


y) (S±b 


y)]F c °f (x»z u ^ 


-z»x v ) 


- 1(Sa 


_A 


2/) (5±b 


x) 


+ (S±A 


x) (S ±B 




-y»z») 


- [(Sj 


_A 


2/) (5± B 


x) 


+ (S±A 


x) (S ±B 


y)]Ffl 2 t>'(x»y» 


+ y tl x v ) 


- [(Sj 


_A 


2/) (5±b 


x) 


- (5_LA 


x) (5xb 






- [(Sa 


_A 


2/) (5±b 


x) 


- (SxA 


x) (Sj_ B 


y)]Fff*'(zV 


+ 2AO 



WET a, 



= (S iA ■ z) (S ±B ■ x) (fi ± + F^') (-<r + 

- 2 (5 M • z) [Sxb ■ x) (F^' + Ff^') x»x v 
+ (S eA ■ z) (5x5 • x) (i& - F e \ - F™ s2 ^') z»z v 

- (S iA ■ z) (S ±B ■ x) F™*'(x»z v + z»x v ) 
+ {StB-y)FtmSt A y»+y»S» tA ) 

+ (5 M • z) Ffl ^ {S%x v + x»SU) ■ (30) 



An analogous relation holds for „, . Notice that, although the complexity that a doubly polarized 

hadron-hadron scattering inevitably presents, at a careful look we can clearly recognize, in each hadronic 
tensor term, the structures already found in the unpolarized and single polarized cases, Eqs. ((U and (fTS]) . 

The general, doubly polarized DY cross section can conveniently be separated in four parts 

da(S A S B ) _ dcr un P da SA da Ss da SASB 
dn<d 4 q ~ dW¥q + dQ,'d i q + dWd 4 q + dWd^q ( ' 
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is the unpolarized cross section of Eq. ([6|), while 



mrsk and are the sin § le Polarized 



where mnF- q 

contributions of Eq. (|17[) . Contracting Eqs. (|28M30j) with the leptonic tensor of Eq. (|A4j) we can compute the 
doubly polarized term of the above cross section 

da SASB a 2 



dQ'd 4 q 2(2tt) 4 M 2 s 
x < S' eA S' eB 



F^(l + cos 2 0') + F«(l - cos 2 0') 

+F™ S *' sin 20' cos 0' + F™ s ^ sin 2 0' cos 20'" 
+S'^ B cost's A + 4f SB ) [Fi x (l + cos 2 0') + F 2 ± (l - cos 2 0') 
+Fl°[ *' sin 29' cos 0' + F?°? 20 ' sin 2 0' cos 20'" 



~S I 4 $ I 



Fj ± (l + cos 2 0') + Fi ± (l - cos 
-Fl°l sin 20' cos 0' + F^ 20 ' sin 2 0' cos 20' 



2 o'\ 



J S A 



J S B 



J S B 



F s ™*' sin 29' sin 0' + fff + sin 2 0' sin 20' 



F^ sin 20' sin, 



y-isin 20' • 2 nl ■ nit 

r ±± sin # sin 20 



S' eA S' ±B cos 0s 



F/ ± (l + cos 2 0') + F 2 j_(l - cos 2 0' 



+F™ S ' sin 20' cos 0' + F™ s 20 ' sin 2 0' cos 20' 
+S't A S'j_ B sin 0g B F™ ^ sin 26' sin 0' + F™ 2 * sin 2 0' sin 20' 



Fl e (l + cos 2 0') + F 2 ^ (1 - cos 2 0') 



-F c ° s *' sin 20' cos 0' + F?T 20 ' sin 2 0' cos 20' 



S'±A S 'eB sin ^ 



Fg 1 sin 20' sin 0' + Ff ? 20 ' sin 2 0' sin 20' 



(32) 



We point out that, as for the case of single polarized Drell-Yan cross section, a transverse/longitudinal spin 
vector in the hadronic cm. frame acquires a longitudinal/transverse component in the dilepton rest frame, 
see Eqs. (|20H22I) . However, as for the case of single polarized Drell-Yan cross section, in the CS frame at 
first order in a qr/M expansion Eq. (|3"2"j) holds with _L= T and £ = L. Eq. (|32p shows a considerably rich 
structure, which is related to all the possible correlations between the hadron and parton spins and the 
parton transverse momenta. 

We can now conclude that the polarized hadronic tensor, Eq. (TT4")) . can be decomposed in terms of 48 
structure functions: 4 for W£% p , 8 for W^, 8 for and 28 for W£" Sb in agreement with Ref. 0. 



IV. PARTON MODEL 



In the previous sections we analyzed the general structure of the Drell-Yan cross section in terms of hadronic 
structure functions. Here, we will explicitly calculate some terms of this cross section in the context of the 
parton model. Ref. Q provides a full classification of the expressions of the (polarized and unpolarized) 
hadronic structure functions in the parton model at born level {qq electromagnetic annihilation), having 
taken into account the intrinsic motion of partons inside the nucleons. 
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Our aim is to show that the decomposition of the hadronic tensor W^ u along the dilepton rest frame 
axes can also be easily applied to the partonic tensor, w l " / , in principle at any order in qr/M (although 
factorization is not guaranteed at any order). At present, we will only consider two relevant examples at first 
order in qr/M: the unpolarized parton distribution function contributions to the unpolarized cross section, 
and the Sivers contributions to the single polarized cross section, which originate from the unpolarized 
partonic tensor. The general expression for the polarized cross section, that requires the full structure of 
polarized partonic tensor, summed over all underlying polarized qq — > £ + £~ annihilations, will be presented 
in a forthcoming paper. 

We work in a kinematical configuration in which hadrons A and B move along the Z crn axis, in the A-B 
cm. frame. Therefore, neglecting hadron masses, their four-momenta are: 



p A = ^(l,0,0,l), p B = ^i(l,0,0,-l), (33) 



where s = (j>a + Pb) 2 - Taking into account the parton transverse motion, the partonic momenta can be 
written as: 

s (~\ 1 k ±b 2 k^ b k ±b 



2 V x b s x b y/s x b sJ 

where x a and x b denote the light-cone momentum fractions of parton a and b respectively. The virtual photon 
7* has a momentum q which, by momentum conservation, equals the sum of the two partonic momenta p a 
and p b : 

q = iPa +Pb) = (<?o,<?t,<7l) with q 2 = M 2 . (36) 

To 0(q T /M) we have: 

„ go + gL go-gi ( r) 

Xa — 1— ) x b — _ \° 1 ) 

V s V s 

which implies x a x b ~ M 2 /s. 

According to the QCD parton model, the hadronic tensor amplitude W^ v is the incoherent sum of the 
corresponding elementary partonic tensors for qq annihilations, w^ v . Therefore, to leading order, the Drell- 
Yan cross section can be written as 



da a 



2 



a^qdVL 1 6M 2 s 



J d 2 k ±a L f q a/SA {x a ,k^ a ) f b/SB (x b ,q T - fc ±a ) , (38) 



where a,b = q,q with q — u,d, s, u, d, s, while L^ v is the leptonic tensor given in Appendix |XJ For the 
time being, we will only consider the contribution of unpolarized partons, therefore the expression of the 
unpolarized partonic tensor in the dilepton rest frame is identical to that of the leptonic tensor, with <f> ! and 
9' replaced by <p' a and 8' a , the azimuthal and polar angles of parton q in the dilepton rest frame 




]j2 + z»z v ) (1 + cos 2 6' a ) - 2 sin 2 0' a z»z v 



+2 sin 2 6' a cos 2<j)' a 



1 fu, q^q" 



2 V M 2 



+ sin 2 6' a sin 2^>' a (x^y l/ + y> l x v ) 



■sm2Q' a cos<i>' a {x^z v + z^x") +sm29' a sm<i>' a {y»z l/ + z^y") } . (39) 
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In Eq. (1381) , f a /s A denotes the distribution function of an unpolarized parton of flavour a inside a hadron A 
with polarization Sa, and is given by the sum of the unpolarized distribution function, f a /A{x a , kj_ a ), and 
the Sivers distribution function, A f a /-f (x a , k a ), which gives the density number of an unpolarized parton 
a inside a transversely polarized hadron A [3, llOj : 

fa/S A ( X a,k± a ) = f a / A (Xa,k±a) + ^A N fa/t( x a,k±a)sm((f>s A -(i)a) 

fb/s B (xb,k± b ) f b / B (xb,k± b ) ■ (40) 
Notice that the Sivers function is also often denoted as fyf(x a , k± a )', this notation is related to ours by 

A N f a/t (x a , k ±a ) = f#(x a , k±a) . (41) 

m p 

The Sivers function is a particularly interesting TMD to unravel the 3-D structure of hadrons, as it embeds 
the coupling among the transverse motion of unpolarized partons and their parent hadron spin, which can 
only be related to their orbital angular momentum. Therefore, a clean extraction of the Sivers function 
will also shed light on the angular momentum of quarks and gluons inside hadrons. Moreover, it is now 



widely accepted that this function should change sign [llj when going from SIDIS to Drell Yan processes, 
a belief that still needs to be supported by experimental evidence. It is therefore of great interest to be 
able to perform a clear and exhaustive analysis of its contribution to Drell Yan processes: although this is 
extremely involved in the general case, as shown in Eqs. (|25[) . it can become particularly simple and manifest 
when the process is examined in the kinematical range in which the photon transverse momentum, qx, is 
much smaller than the dilepton invariant mass, M , and of the same order of magnitude of the quark and 
anti-quark intrinsic momenta, as we will illustrate in what follows. 
Contracting the unpolarized leptonic and partonic tensors we obtain: 

L^w^ = — |(l + cos 2 6>')(l + cos 2 6K) + 2 sin 2 6»'sin 2 ^ 

+4 sin 2 9' cos 20' sin 2 9' a cos 20' o 
+ sin 2 9' sin 20' sin 2 0' a sin 2cj)' a 
+ sin 29' cos 0' sin 29' a cos (j)' a 

+ sin 29' sin 0' sin 20' a sin <j>' a \ . (42) 

Inverting Eq. (|B4|) we can connect the components of partonic momenta in the dilepton rest frame with 
those in the hadronic cm. frame; to first in order qx/M we find: 

cos 2 ^ ~ 1 (43) 

sin 2 6»l ~ (44) 
2 

sin 26>^ cos </> a ~ - — (q T - 2k ±al cos0 7 - 2k ±a2 sin0 7 ) (45) 
ivl 
4 

sin 29' a sin 4> a ~ - — (fc X ai sin 7 - k ±a2 cos cf> 7 ) , (46) 

where k± a \ and k± a2 are the transverse components of p a in the hadronic cm. frame. Eqs. (|43[) and (|44l) 
tell us that only three leading terms survive in Eq. (|42l) . Thus: 



L""uv ~ MH (l + cos 2 e') 



, / It 2fc ±a i 2k^ a 



sin 29 cos <p \ cos <p^ sin 

MM M 



sin 29' sin 4>' ( " , "~ Lal sin 7 cos 7 ) ^ . (47) 

M M 
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At this stage, the final expression for the cross section is obtained by inserting Eqs. (|4"0)) and (14"T|) into 
Eq. (|38|) . The integration on k± a can be further simplified by reconsidering the integrals after a rotation of 
the system to the "virtual photon" production plane, as described in detail in Appendix [Cl Finally we find: 



da 



d 4 qdfl> 6M 2 s 



cos 2 9' 



- sin 29' cos cj>' d 2 k± a 



2 — (k^ a -q T )-- 



+ (1 + cos 2 9') sin(4>s A - 4> 7 ) I d z k ±a (k ±a ■ q T )A N f« n (x a , k ±a )f« /B (x b , k Xb ) 



■sm20'cos0'sin(<^ - 7 ) / d k± a (kj_ a ■ q T ) 



M 



+ sin 29' sin 6' cos( 



) / d'k ±a 2- 



A I 



(fe±a • 9t) 2 - l| A N f g a/t (x a ,k± a )f^ /B (x b ,k ±b ) ) , (48) 



where we recall that the lepton angles 9' and </>' are expressed in the CS frame, while 4>s A 
to the hadronic cm. frame. This relation deserves a careful discussion. 



and 6-y are referred 



• In the first line we find the leading term of the unpolarized cross section, proportional to the convolution 
of the two unpolarized parton distribution functions, f a /A{ x ai k± a ) ® fb/B{ x bi k±b), and identified by 
the (1 + cos 2 9') angular distribution. 

• In the second line, we find the first order correction to the previous term in our qx/M expansion 
(notice that k±/M and qr/M are of the same order). This term can be considered as the analogous 



of the so called Cahn effect in SIDIS 12|, |l3( . In Eq. (IC9[) of Appendix [C] we show that it is directly 



proportional to the difference between the mean average transverse momenta of parton a and of parton 
b inside the two colliding hadrons. In the Gaussian model hypothesis, see Eq. (|52|) . this term translates 
directly in the difference between the Gaussian widths. Thus such term is negligible when the two 
colliding hadrons are charge conjugated, as it is the case for pp induced DY processes, while it could be 
observed in np DY processes if the pion and proton transverse momentum distribution are different, 
or in pp scattering if the quark transverse momentum distribution (valence contribution) is different 
from that of the anti-quark (sea contribution) . Although the existence of this term in the unpolarized 
cross section was pointed out many years ago in Refs. [l| and [7|, it is only now, in light of the 
new theoretical developments and of the planned experimental programs, that its phenomenological 
implications related to its detailed form are becoming particularly relevant. Finally, we point out 
that here we do not find the term proportional to sin 2 9' cos 20', i.e. the so-called Boer-Mulders term, 
which is usually included in the unpolarized DY cross section. This is due to the fact that we are not 
considering the contributions originating from any polarized elementary qq annihilation. Its absence, 
however, does not invalidate our phenomenological considerations, as this term gives rise to a different 
azimuthal dependence and can, therefore, be experimentally isolated. 

The last three lines correspond to the contributions of polarized hadrons to the cross section, and involve 
the convolution of the Sivers distribution function A N f a /f(x a , k± a ) with the unpolarized distribution 
function fb/B(xb, k± b ): line three corresponds to the dominant Sivers effect which, in DY processes, is of 
great importance from the theoretical point of view. As mentioned above, the Sivers function reverses 
its sign when is observed in DY and SIDIS processes l_l| : presently, no polarized DY experiments are 
able to confirm this prediction but many dedicated future experiments are planned in several facilities: 
RHIC at BNL 0, PANDA [H] and PAX at FAIR, COMPASS at CERN IHEP [3] in 
Protvino, and JINR 19] in Dubna. 
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• The last two lines of Eq. 
... sin 26' wa.{(j)s A — 



, which can be rewritten in an alternative, more suitable form as 
1 



+ sin 26' sm(<ps A - <£ 7 - <!>') 



d z k. 



d 2 k. 



2M 
2*i. 



4fc ±a (fe ia • q T f 
- q T (k± a ■ q T ) 



k±a ■ (q T - 2fe_L a ) 



Ja/f J b/B 



2M 



A N f q f q 

^ Ja/t Jb/B > 



(49) 



represent the equivalent of the "Cahn effect" for the single polarized DY, and show a structure very 
similar to that of the unpolarized part. These contributions are suppressed by one power of qx/M 
and, as in the unpolarized Cahn effect, we can recast the integration on k± a similarly to what is done 
in Eq. (|C9[) : by doing this one can observe that they can give access to the difference between the 
average transverse momenta of the unpolarized and the Sivers distribution functions. 

These results become particularly transparent and acquire a phenomenological value when the integrals over 
the intrinsic transverse momenta are explicitly performed by using a simple Gaussian model for the TMDs. 
Similarly to what was done in the last Section of Ref. [l(| , we assume the k± dependence of the TMDs can 
be factored and approximated with a Gaussian distribution of the form: 



fa/A(x a ,kj_ a ) = f a / A (x a ) 



(50) 



where f a /A ( x a ) , can be taken from the available fits of the world data. In general, we allow for different widths 
of the Gaussians for the different parton flavours, but take them to be constant. For the Sivers function, we 
assume a similar parametrization, with an extra multiplicative factor k± a to give it the appropriate behavior 
in the small k± a region [2oj |: 



Af a/t (x a ,k± a ) = A/ Q/t (x) v2e — ^ — 



(51) 



where the x-dependent function Af a /^(x a ) is not known, and should be determined phenomenologically by 
fitting the available data on azimuthal asymmetries and moments; the k± dependent Gaussian has been 
assigned a width ( fc 2 ) s and a suitable normalization coefficient \[2e to make sure it fulfills the appropriate 
positivity bounds [21] . 

By inserting Eqs. ((50]) and ([5TJ into Eq. ((48]) we get 



d(J unp 



ST J 



d 4 qdfl' 6M 2 s ^ €q r °/^ Xa) f "/B 



e -«r/<«r) r 
( x b) _ / „ 2x (( 



t<9t) 



1 + cos 2 6') 



and 



dcr SA 



2ea 2 



d 4 qdW 12M 2 s M s 



— ^4 ^f q an M fl /B (x b ) 



M ( 



-9t/(9t>< 



(k 2 



It) 



• sin 26' cos 



(52) 



<*5J 



| (1 + cos 6') sin(05 A - 



qr_ 

' M 



q\ 



(klb) 



(l T )s 
sin 26' sin 



(k±b) 



sin 26' cos d> sin( 



<i>' cos(4> Sa - <t>l) } 



where we have defined (kj_ a 
rearranged as Eq. ((25]) : 



> + (kl b ) = (g 2 ), and (kl a )l + (kl b ) = (q 2 ), 



(53) 

Notice that Eq. ((53)) can be 



da SA 
d i qdVL' 



y/2ea 2 qr sr^ 
12M 2 s ~M~S^ 



e\ Af a/t ( Xa ) f b/B (x b ) 



b/B\ 



n((q? 



Tl s i 
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|(1 + cos 2 6')sm((f)s A - </> 7 ) 
+ sin 20' sin((f> SA — </> 7 + </>') 



-— sin 29' sin^s^ — </> 7 



q T (klJ a - (fei 6 ) 



2M (&) g 



(4) 



Eqs. ([53")) and (|54p have a great phenomenological value as they show that if a sin(^s A — <f> 1 
modulation is experimentally detected, then it gives direct access to the asymmetry (fc^_ a ) s - 



(k 



(54) 



b') azimuthal 
5_ 6 ), whereas 



the sin0 cos(4>s A 



azimuthal modulation would give access to (kj_ b ) alone, and would survive even in 



the case in which (kj_ a ) s — (k\ b ). 

Notice that, from SIDIS experiments we know that (kj_) s ^ in particular, in the latest data 

fits 2(J 22| the ratio (k±) s /(k'j_) is rather well constrained and turns out to be between one half and two 
thirds. Consequently, the scenario of equal fcj_ distributions for the unpolarized and Sivers distribution 
functions seems to be excluded, making even more interesting any further investigations in DY processes, 
where we could have access to the values of (kj_) s and (kj_) separately. 



CONCLUSIONS 



The theoretical and phenomenological relevance that Drell-Yan processes have recently assumed in the 
context of hadron physics is out of any doubt. New dedicated experiments are either under construction 
or beingplanned, like RHIC at BNL PANDA [H] and PAX [16] at FAIR, COMPASS at CERN 
IHEP 18] in Protvino, and JINR [l9} in Dubna. They will help to shed light on many issues, in particular 
on our understanding of the 3-D dynamical structure of nuclei and, more in general, of hadrons in terms of 
their elementary constituents, quarks and gluons. This motivates our efforts toward the full understanding 
of the DY general cross sections, understanding that is mainly achieved by the decomposition of the hadronic 
tensor in terms of structure functions. 

We have performed this decomposition by using the helicity structure functions, i.e. projecting the hadronic 
tensor along the dilepton helicity axes. Some work was already done in this direction 0,SI13] : m particular, 
m Ref. the hadronic tensor was decomposed with a tensorial basis constructed with the dilepton momen- 
tum, the hadronic momenta and the spin vectors, which resulted in a decomposition in terms of invariant 
structure functions. Although this basis leads to the desired results, it is not ideal to perform such a decom- 
position, as it generates a huge number of non-independent structure functions that must subsequently be 
eliminated by exploiting all the symmetry properties the hadronic tensor should exhibit. 

In our work we decompose the hadronic tensor using the dilepton helicity axis as tensor basis, a choice which 
ensures that all the symmetry properties are fulfilled from the beginning, avoiding the initial proliferation 
of non-independent structure functions. Moreover, the helicity structure functions can be intimately related 
to the elementary scattering process, being connected to the polarization status of the virtual photon. 
Therefore, we believe we have found a suitable and elegant working setup which enables us to disentangle 
the 48 structure functions of the polarized DY cross section in the most general way, in a configuration valid 
in any dilepton rest frame, with a minimal amount of effort. 

Finally, within the QCD parton model approximation, we have performed a phenomenological analysis of 
the unpolarized term of the DY cross section and of the Sivers effect, both originating from the elementary 
scattering of an unpolarized qq pair, in the kinematical range in which kj_ ~ qx <C M. Interesting results, 
which could help gathering valuable information on the k± distribution of the unpolarized TMD distribution 
function, f q / p (x q , k±), and of the Sivers TMD distribution function, A N f q / p ^ (x q , k±), are presented 



The qr/M suppressed terms in Eq. ([48)) deserve a similar phenomenological analysis, that we are planning 
to perform in a forthcoming paper. One can also wonder if they could be studied above the so-called safe 
region, where the suppression should be less pronounced: in principle, if in a single polarized DY process 
the J/ip is found to be transversely polarized (in the CS frame) one can apply Eq. (|4"5)) . with obvious 
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modifications to accommodate for the J/ip coupling, in analogy to what has been done in Ref. 24|. A 
further extension of this work, including the evaluation of the terms generated by the transverse polarization 
of partons, therefore related to the Boer-Mulders and transversity functions, will soon be presented. 
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Appendix A: Leptonic Tensor 



For reference we report the structure of the unpolarized leptonic tensor decomposed by means of the 
helicity axes, see Ref. [23[. It can be easily obtained by inserting the expressions of 1% and l\ in terms of the 
hclicity axes into Eq. ((2]) 



/'' 



/'' 



{qf* + M (sin 6' cos (j)'x^ + sin 6' sin 4>'y^ + cos O'z*)] 
[cf - M (sin 9' cos ftx^ + sin & sin V + cos & z»)\ 



and keeping in mind that 



as x M , y M and z^ are all perpendicular to q^. Thus one obtains: 

M 2 ' 



L» v = — 



r 



q^q v 



2 sin 2 9' cos 20' 



V z"z v \ (1 + cos 2 9') - 2 sin 2 9'z^z" 



q^q v 



+ z fi z" 



+ sin 2 9' sin 2<j)'(x ti y v + y^x") + sin 29' cos <j)'{x^z v + z^x v ) 
+ sin26>'sin0'(?/ 1 z ly + z'V) \ ■ 



(Al) 
(A2) 

(A3) 



(A4) 



Appendix B: Collins-Soper helicity frame 

The leptonic angular distribution is studied in the dilepton rest frame where leptons are back to back. 
The dilepton rest frame can be chosen in many different ways. A very common choice is the so called 
Collins-Soper (CS) frame. 0- In the hadronic cm. frame the CS axes read: 

z CS - . ' (gL.O.O.ft) (Bl) 



VM 2 + It ' 



— — , COS U>~ , Sill w~ , 

y/AP + q z T V M ' M Vl1 M Vl1 M ) 

(o, — sin0 7 , cos0 7 , 0^ (B3) 



15 



where = (go, Qt cos </> 7 , qr sin 7 , q^) is the dilepton (virtual photon) momentum in the hadronic cm. 
frame. Notice that the y cs axis has no time-like nor longitudinal components. 
The matrix: 



/ .22. 
/ M 



AM" 



(JT COS (f)-y 

M 
qT sin 



M C0S 

1 V^+g^ 



V 



M 
A/ 



M 



_2M2_ 







COS(^ 





1L 









go 



(B4) 



is the Lorentz transformation from the Collins-Soper frame to the hadronic cm. frame. Notice that the 
columns of this matrix are the four column vectors f = q^/M, , y M , in the Collins-Soper frame. 



Appendix C: Integration by rotation in the virtual photon production plane 

Let us define as the "virtual photon" production plane the plane containing the virtual photon momentum 
q and the Z cm axis in the hadron cm. frame. We can then define a new frame where = Z cm and the 
X-y axes lay in the virtual photon production plane and Y 7 is perpendicular to both X-y and Z 7 . This new 
frame is rotated by an angle 7 with respect to the cm. frame (X cm , Y cm , Z cm ) : 

X Cm = Xry COS ^ *y Sin (j)ry (CI) 

^ Cm = Xfy Sin (f)j + V-y COS (fi-y (C2) 

Z cm = ^7 (C3) 

As a direct consequence, the transverse photon momentum will be along X 7 = q T . Notice also that, 
for partonic momentum conservation, any function of the parton intrinsic transverse momentum moduli, 
9(k\ al k\ b ), obeys the relation g{k 2 j_ a , k 2 j_ b ) = g[k ±a ,(q T - k ±a ) 2 } = g(k 2 La ,k ±a ■ q T ). When we insert 
Eq. (|47|) into Eq. f|38[) the first two terms, proportional to the convolution of the unpolarized distribution 
functions f a /A fb/Bj will show integrals of the form: 

J d 2 k ±a g(k 2 La , k ±a ■ q T ){k ±al cos0 7 + k ±a2 sin0 7 ) 
= d 2 k ±a g(k ±a ,k Xa - q T )[(k_L 

a ' X cm) 

cos (f)ry + (k± a ■ Y cm ) sin c^ 7 ] 
= / ' d 2 k ±a g{k 2 ±a ,k ±a - q T ){k ±a - X 7 ) , (C4) 

/d 2 fe_L a p(fci a ,fe_L • q T )(k± al sin 7 - kj_ a2 cos^ 7 ) = - / d 2 fc ±a g(k 2 ±a ,kj_ a • q T )( k -i-a • Y" 7 ) = 0, 

(C5) 

having exploited the fact that J d 2 k± a g(k 2 _ a , k± a ■ q T )(k± a ■ V 7 ) = 0, because (fej^ a • Y" 7 ) is an odd function 
over the symmetric range of fcxa integration. 

The Sivers terms are more involved, as the Sivers function itself depends on kj_ a through its proper phase, 
sm((/)s A — 4>±a) oc (sm(j)s A k± a i — cos(f>s A k± a 2)- Therefore, we will have to deal with integrals of the form 

J d 2 k ±a g(k 2 ±a ,k ±a ■ q T )(k± a i sin0s A - k ±a2 cos(/) Sa ) 

= J d 2 k ±a g(k 2 La ,k ±a ■ q T )[{k ±a ■ X 7 ) sin(0g 4 - 7 ) - (k ±a ■ Y 1 )cos((f>s A - 7 )] 

d 2 k± a g(kj_ a , k ±a ■ q T )(kj_ a ■ X~) sin(0 5A - </> 7 ) , (C6) 
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J d 2 k^ a g(k 2 ±a ,kj_ a ■ q T )(kx a i cos0 7 + kj_ a2 sin0 7 )(fcj_ Q i sin^ ~ fc ±a2 cos^s A ) 

= J d 2 k ±a g(k 2 Xa , k ±a • q T )(fe_L Q • X 7 ) 2 sin(0 SA - 7 ) , (C7) 

J d 2 k^ a g{k 2 _ al k^ a ■ q T )(kj_ al sin<^ 7 - k^ a2 cos c/>^)(k± al sm^s A ~ k ±a2 cosc/>s A ) 

= J d 2 k ±a g(k 2 ±a ,k ±a ■ q T )[k\ a - (k ±a • X 7 ) 2 ] cos(0 Sa - <£ 7 ) . (C8) 

Finally, it is very interesting to notice that 

J d 2 k± a g(k 2 ±a , k± a ■ q T )[qT - 2(k± a ■ X 7 )] 

= / d 2 k ±a d 2 k ±b S 2 (kj_ a + k ±q - q T )g(k 2 ±a ,k 2 ±b ) — [(k ±b - k ±a ) ■ (k ±b + k ±a )} 
J It 

= f d 2 k ±a d 2 k ±b S 2 (k^ a + k u - q T )g(kl a , fc 2 6 ) fc lb- fc ia ; (C9) 



which lead to the conclusions outlined in Sect. [TV 
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